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Abstract: We compute the fully-differential B —* X U W decay width to all orders in 
perturbation theory in the large-/?o limit. Each of the five structure functions that build 
the hadronic tensor is expressed as a Borel integral, summing up 0(Cf/3q _1 a™ ) corrections 
for any n. We derive analytic expressions for the Borel transforms of both real and virtual 
diagrams with a single dressed gluon, and perform an all-order infrared subtraction, where 
the Borel parameter serves also as an infrared regulator. Expanding the result we recover 
the known triple-differential NLO coefficient, and obtain an explicit expression for the 
0(CfPoo?) triple-differential NNLO correction. This result can be used to improve the 
determination of from inclusive B — > X U W measurements at the B factories with a 
variety of kinematic cuts. 
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1. Introduction 

The measurements [1,2] of semileptonic b decay branching fractions (BF) play a crucial 
role in the determination of the CKM matrix elements [3,4], which form the basis for many 
precision tests of the Standard Model and provide an input for new physics searches. While 
any potential discovery of new physics in the flavor sector is associated with loop-induced 
transitions, the CKM parameters are most reliably determined by tree-level weak decays. 
Here two fundamental ingredients are |V^b| and |Vub|, which are measured in semileptonic 
b — > c and b — > u decays, respectively. 
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Both inclusive and exclusive semileptonic measurements are used to extract these pa- 
rameters. Inclusive measurements are inherently more robust owing to their limited sen- 
sitivity to the hadronic structure of the initial and final states. However, since b — ► u 
transitions are about 50 times less abundant than b — > c ones, kinematic cuts must be ap- 
plied in order to isolate the b — ► u decays and measure | V^b | • Consequently, the calculation 
of the fully differential spectrum is essential for precision measurements of | V^b I • 

The theoretical calculation of inclusive decay spectra is complicated by the presence 
of large perturbative and non-perturbative corrections. In b decays into light quarks, 
e.g. B — > Xgj and B — > X U W, most events are characterized by jet-like momentum 
configurations, where the invariant mass of the hadronic system in the final state is small. 
When computing the differential spectrum, or the BF with kinematic cuts, one encounters 
parametrically-large Sudakov logarithms as well as non-perturbative corrections that are 
associated with the momentum distribution of the b quark in the meson [5-14]. 

Recently, there has been significant progress in the application of resummed pertur- 
bation theory to compute inclusive decay spectra using the method of Dressed Gluon 
Exponentiation (DGE) [15-19]. Underlying this approach is the realization that running- 
coupling corrections play an important role in shaping the spectrum. Beyond the purely 
perturbative issue, infrared renormalons are useful in consistently separating between per- 
turbative and non-perturbative corrections while retaining the predictive power of pertur- 
bation theory. 

The significance of running-coupling corrections stems from the fact that the gluon 
virtuality, which sets the effective scale of the coupling [20], is typically lower, sometimes 
significantly lower, than the hard scale m^ that is used as the default renormalization point. 
Consider for example the fully-integrated b — ► X u lv width, 

\ 2 



X u lu)= G2 ^ Vuh?m " h 



192vr 3 



which is known to NNLO [21], 
b\ = C F [ tt z 



1 + h Os ^ +b J^ W + 
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Po = ^C A -±N f , (1.3) 

where Nf is the number of light flavors. Here we have split the 62 coefficient computed in 
Ref. [21] into a running-coupling part, proportional to (3q, and the rest. We find that with 
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Nf = 4 the former yields b^ — —26.84, while the latter &2 6St — 5.54, adding up to 62 — 
— 21.30. Evidently, the running-coupling corrections are dominant. These large corrections 
are related to the leading infrared renormalon of the pole mass nib, which is located at 
u = 1/2, where u is the relevant Borel parameter. The eventual cancellation [22-24] 
of this ambiguity in Eq. (1.1) involves the overall factor m\ on the one hand and the 
series in the square brackets on the other. This means, in particular, that higher-order 
corrections 0(C,p/3Q _1 a") in Eq. (1.1) are large and form a non-summable series. Owing 
to the proximity of the u = 1/2 renormalon to the origin of the Borel plane, and the 
relatively low scales involved, the factorial divergence becomes relevant already at the first 
few orders. 

Let us consider now the differential decay width. As usual, the effective scale de- 
pends on the kinematics; in the region selected by kinematic cuts, where the invariant 
mass of the hadronic system is small, radiation is confined to be soft or collinear. An 
obvious consequence is that the effective scale of the coupling gets small, and therefore 
large running-coupling corrections should be expected. 

Moreover, the normalized spectrum too is affected by infrared renormalons. Despite 
the absence of the overall factor m|, infrared renormalons show up in the normalized 
spectrum because the support of the on-shell decay width is set by m&: an 0(A) variation of 
the pole mass amounts to an 0(A) shift of the Born-level (5-function spectrum. Therefore, 
all the moments of the normalized b decay spectra, defined at the partonic level, have 
an infrared renormalon ambiguity at u = 1/2 [25]. In DGE this ambiguity, which affects 
the Sudakov factor, gets cancelled against the pole mass upon computing the resummed 
spectra in physical, hadronic variables [16]. In other approaches [26,27] it is dealt with by 
using an infrared cutoff and absorbing the soft contribution into the definition of the non- 
perturbative parameters. In any case, the presence of this infrared sensitivity at the level 
of the partonic calculation cannot be ignored. Computing decay spectra to higher orders 
in perturbation theory, one therefore expects to find large running-coupling corrections. 

Recently, a first complete NNLO calculation of an inclusive decay spectrum has been 
performed [28,29] for the case of B — > X s ~f through the effective magnetic dipole operator. 
A striking feature of this result is the dominance of the 0(CfPoo^) contribution (which 
has been known since a while [30]) with respect to other color factors appearing at this 
order. The similarity of the two processes, B — > X s j and B — ► X U W, and the dominance 
of running-coupling corrections in the former, leave no doubt that these corrections are 
dominant also in the latter. 

In the case of the triple-differential B — ► X u lu spectrum, the perturbative expansion 
is known in full to 0(a s ) (NLO) only [31]. NNLO corrections have been computed in full 
only for the integrated width [21], Eq. (1.1) above. In addition, 0(CfPoo^) real-emission 
corrections have been recently computed [32] for one particular single-differential spec- 
trum, namely the distribution in the (small) "plus" lightcone-momentum component, pf. 
The complete 0(Cf(3qOi™ +1 ) have also been computed numerically for the five structure 
functions of B — ► X c lv [33]; the results for physical observables obtained with a finite 
charm mass can in principle be extrapolated numerically to the massless case, but this 
procedure involves delicate numerical cancellations and lacks the flexibility necessary in 
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practical implementations. 

Additional results beyond the NLO are available in the Sudakov limit [15, 16]: the 
Sudakov exponent has been determined at NNLL accuracy [34-36], and to all-orders in 
the large-/3o limit [16,25]. 

In this paper we perform an all-order calculation of the triple-differential B — > X u lv 
spectrum in the large-/3o limit. We derive analytic expressions for the Borel transform 
of real and virtual diagrams with a single dressed gluon, which represent the sum of 
0(Cf(3q corrections for any n. We then preform an all-order infrared subtraction 

directly in terms of the Borel variable. By expanding the result we recover the known 
triple-differential NLO coefficient [31], and obtain an explicit expression for the 0(Cf(3qoi 2 s ) 
triple-differential NNLO correction. By integrating this expression we confirm the results 
of the single-differential p~j~ spectrum [32] as well as the (3$ term in the integrated width [21]. 

The 0(Ci^/3oCKs) triple differential width we compute here is an important ingredient 
in improving the determination of [V^bl from inclusive B — > X u lu measurements with a 
variety of kinematic cuts. In this paper we do not perform any numerical studies; these 
will be reported on separately. 

The structure of the paper is as follows. In Sec. 2 we recall the kinematics and set up 
the notation. Next, in Sec. 3 we present the real-emission "characteristic function" based 
on the calculation of Ref. [33], which was performed with an off-shell gluon. In Sec. 4 we 
derive a Borel representation of the real-emission corrections; in Sec. 5 we expand the Borel 
function to obtain explicit formulae for coefficients at the first few orders. Next, in Sec. 6 
we consider the Sudakov limit and extract the non-integrable terms in the real-emission 
contribution to all orders, confirming the results of [25]. Using the Borel variable as an 
infrared regulator, we prepare the tools for an all-order infrared subtraction. In Sec. 7 we 
compute the virtual diagrams, using a Borel-modified gluon propagator. We then perform 
the subtraction of infrared singularities, directly in terms of the Borel variable. In Sec. 8 we 
combine the results of both real and virtual diagrams for the different structure functions, 
getting explicit expressions for the coefficients of the triple-differential rate at NNLO. In 
Sec. 9 we demonstrate the way in which infrared subtraction gets modified depending on 
the kinematic variables used. Finally, in Sec. 10 we summarize our conclusions. 

2. Definitions and kinematics 

Let us write the triple-differential width of 

b(p) —> l(pi) + P(pv) + X u (pj) 

as 

_dr = G^IVubl P )wHp, q), (2.1) 

dpjdpjdEi 16vr 3 w,<ih V ) 

where the total momentum carried by the leptons is q = p\ + pp and E\ is the energy of the 
charged lepton in the b rest frame. The momentum of the hadronic system is expressed in 
terms of lightcone components, namely 

pf = E jT \pj\, (2.2) 
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so p- = ami, and p^ = f3m,f } are the large and small lightcone components of the "jet", 
respectively. They obey 

O<0<a<l. (2.3) 

The relation of these variables with the invariant masses of the "hadronic" (partonic) and 
the leptonic systems, respectively, is 

p) = m\afJ- q 2 = m 2 (l-a){\-(3), (2.4) 

where < p 2 < ^1 — \f<f^j ■ Throughout the paper, represents the pole mass of the b 
quark. 

It is convenient to write the normalization of the differential width in terms of the 
total tree-level width 



namely express Eq. (2.1) as: 



d 3 T 1 dT 1 

, ,„ , = - + - _ — = 6 To — L^( Pl ,p p )W^(p, q), (2.6) 
dadfidx 2 dpj dpj dE t m£ 

where x = 2Ei/mi J . Phase-space integration then yields: 

rl rl rl-ra 

Fb^u = da dra / dx 

Jo Jo Ji-a dadfidx 

rl rl rl—ra ^ 

= 6T da dra dx — L^p^p^W^ (p, q), (2.7) 

Jo Jo Jl-a m b 

where we use r = (3/a following [16]. 

The choice of the lightcone variables in Eq. (2.1) is motivated by the fact that the 
final state is typically jet-like: at Born level p 2 = so most events are characterized by 
p 2 <^m 2 , namely (?«a. While the first NLO calculation of the triple-differential b — > ulv 
spectrum [31] used other kinematic variables, the advantages of lightcone variables have 
recently been acknowledged by several authors [13,16,26,32]. 

For massless leptons, the leptonic tensor is given by: 

L^(pi, P p) = P fp u p + pH ~ Pi ■ PugT ~ i^ pv °Pl P Pv*- (2-8) 

Lorentz decomposition of the hadronic tensor W /MI/ (p,q) gives rise to five scalar "structure 
functions" : 

W">,9) = -Wi(a,j9)r + W 2 (a,p)v> i v v + iW 3 (a, f3)e fiupa v p q a 
+ W 4 (a,(3)q^q u + W 5 (a,/3)(^<f + v"$ i ), 

where v = pfm^ and q = q/rrib- W J,u (p,q) is related to W^ u (p, q) defined in [33] by: 



d(q 2 , P 2 ) 



d(pj, p- 



(« Pl W ^(p,q). (2.10) 



7T 
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Note that both Wf v (a, (3) and Wf v (a, (3) are dimensionless. Contracting the Lorentz 
indices between the leptonic and hadronic tensors, Eqs. (2.8) and (2.9) respectively, one 
finds: 



l -L^( VuPu )W^( V ,q) = (l-a)(l-/J)Wi(a,/?) 



- ^{x 2 - x (2 - a - (3) + (1 - a)(l - 0)j W 2 (a, /?) (2.11) 

+ (1 - a)(l - 0) (x - 1 + \{a + /?)) W 3 (a, (3), 

where, as above, x = 2Eijm^. Each structure function Wj(a,/3) may be decomposed as: 

Wi(a, /?) = ^(a) <S(/3) + R^a, (3). (2.12) 

The functions Vi(a) and Ri(a,f3) have perturbative expansions in a s (m^). At the leading- 
order (LO), Ri(a,/3) = and, for i = 1 to 5, 

^ LO (a) = [a, 4, 2, 0, -2] . (2.13) 

Substituting (2.12) with (2.13) into (2.11) and using (2.6) one gets: 
d 3 T LO 

——— = T uj (a,x)5((3); u (a, x) = 12 {x + a - 1) (2 - x - a) (2.14) 

da dp dx 

where To is given in (2.5). As usual, this Born-level result receives perturbative cor- 
rections from both virtual and real-emission diagrams. Purely virtual contributions are 
proportional to 5(f3). These, however, contain infrared singularities that cancel against soft 
and collinear real-emission singularities when performing phase-space integration near the 
[3 — ► limit. Thus, beyond the LO, the separate definition of V$(a) and Ri(a, (3) in (2.12) 
requires a subtraction prescription. We shall return to this issue in Sec. 6. 

The full NLO, 0(a s ) result has been obtained in Ref. [31], and checked in various pa- 
pers that considered higher-order running-coupling effects, including Ref. [25] (real emis- 
sion) and Ref. [33] (real and virtual corrections). The result was presented in terms of 
lightcone variables in Ref. [16]. 



3. Real emission of an off-shell gluon and the characteristic function 

Perturbative calculations of many observables in QCD, including inclusive cross sections 
and decay rates can be improved by the resummation of running-coupling effects [20, 33, 
37-42]. Specifically, keeping just the leading term in the (3 function, one sums up the terms 
proportional to /3Q _1 a™ to all orders, the so called BLM terms. 

Technically, running-coupling terms can be conveniently derived using the dispersive 
method, see e.g. [38,41,42], where the one-loop calculation is performed with an off-shell 
gluon. The calculation of the semileptonic decay "structure functions" Ri(a,f3) with a 
single off-shell gluon was performed in Ref. [33] for the more general case of b — > c decay, 
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where the charm mass m c is kept. Here we use this result to derive the corresponding b — > u 
result by sending m c — ► 0. This limit leads to significant simplification that facilitates 
obtaining closed form analytic expressions. 

Further simplification is achieved by choosing the lightcone variables described above, 
which are suitable for the final state "jet" kinematics. As we shall see, the result for Ri(a, (3) 
is completely symmetric under a < — ► (3; only the phase-space restriction, Eq. (2.3), breaks 
this symmetry. 

The LO calculation of the real-emission diagrams with a gluon of virtuality m 2 = £,m 2 
yields: 

R t (a,P) ^C F ^™^F t (a,P,0 + O(a 2 s ) (3.1) 

TT 



with the following "characteristic function" : 



(a(3) Zi \t- t+ J iy " a/3 \t 



(3.2) 



where the powers are yi = [1, 2, 1, 2, 2] and = [3, 3, 3, 2, 3] for i = 1 through 5, respectively, 
and Pi(a,/3,£) and Qi(a,(3,£) are polynomials in all their arguments. Finally, 

where 

v 7 ! = a - (3. (3.4) 
Note that the phase-space limits are 

< £ < a(3, (3.5) 

where the upper limit corresponds to the situation where the entire mass of the hadronic 
system p 2 - = m 2 a(3 is given by the gluon virtuality. 

4. Borel representation 

The Borel representation of the result for Ri in the single-dressed-gluon (SDG) approxi- 
mation is: 

RT G {a^) = duT{u) Br G (a,(3,u) (4.1) 

Po Jo \ m bJ 

n \ Wt n\ a s( m b) . (2), a\a f a s {m b )\ 2 
= C F c\'{a,[3) he- J (a,/3)(3 [ +■■ 

7T \ TT J 

where mi, is the bottom pole mass, A and a s are defined in the MS scheme, and (3$ is defined 
in Eq. (1.3). T(u) is the inverse Laplace transform of the coupling (see Appendix A); in the 
large-/?o limit, where the renormalization-group equation is just one loop, T(u) = 1. Re- 
summation of running-coupling corrections beyond this strict limit can also be performed 
using (4.1). This is briefly explained in Appendix A. 
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The Borel function in Eq. (4.1) can be derived from the following integral over the 
characteristic function (see e.g. [41,43] or Sec. 2.2. in [44]), 



B?° G (a,(3,u) ee - e §«™ f ^(a, (3, 

^ Jo Q 

= _jus™™ ( a( 3y u bi (a,(3,u), (4.2) 

7T 

where jFj(a,/3, £) is given in Eq. (3.2) and, upon changing variables from £ to r, where 
£ ee (1 - r,)a(3, 

bi(a, (3, u) ee /' dr, (1 - r/)- 1 -" ^(a, /?, (1 - r?)a/3). (4.3) 
J o 

We have performed the integral in (4.3) analytically, and checked the result by numerical 
evaluation. Below we give a few details of the calculation and summarize the analytic 
expressions for bi(a, (3,u). 

Writing of Eq. (3.3) in terms of rj we have = (a(3 — a — (3)(1 — K±rf), with 

±_ a(3 -a- (i [ ' 

so using Eq. (3.4) n + = (a/3 — a) /(a/3 — a — (3) and k_ = (a(3 — f3)/(a(3 — a — (3). Given 
Eq. (3.2), the basic integral needed in Eq. (4.3) is of the form 

f d * - ">~"r^ = ^ ^. (Ml, P - »], « ± ). («> 

All the terms in Eq. (4.3) can be expressed in terms of this specific hypergeometric function 
with the two assignments of the argument, k±, and some additional rational functions. For 
example, to integrate the log term in Eq. (3.2) times (1 — r,) 3 where j is a positive integer 
(to account for Qi(a,(3,£) that are quadratic in £ we need j = 0, 1,2) one first integrates 
by parts and then uses Eq. (4.5) to obtain: 



I dr, (1 - r,)- 1 -^ ln(l - k ± t]) = — — [ dr, (1 - r,) j 
Jo 3 ~u Jo 1 - K ±V 



-K± 



(J ~ u)U ~ u+ 1) 



- 2 Fi([l,l],[2-u + j],/e ± ). 



(4.6) 



Finally, there are known identities that facilitate integer shifts of the indices of hyperge- 
ometric functions. For example, to express the hypergeometric function in Eq. (4.6) in 
terms of our basic one in Eq. (4.5) we need to shift the third index from [2 — u + j] into 
[2 — u\. This is straightforward to do using Eq. (2.10) in Ref. [45]. 
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The final result for bi(a,f3,u) takes the form 



bi(a, (3,u) = A 



Di(a, (3, u) ( 2 F 1 ([1, 1], [2 - u], k+) - 2 F X ([1, 1], [2 - u], 
+ VA$(a,/3,u)( 2^1 ([1,1], [2 -«],«+) + 2 Fi ([1,1], [2 



+ y/\Ti(a,P,u) 



(4.7) 



where the entire u dependence of the coefficient functions is summarized by 

, , A,i(«,/?) , DiMJ) 

U ^P^ U )- u + !_„ + (1 _ n)2 + 2 _ n 

c , fl x ^,0(0,/?) , Sn(a,(3) Sn(a,(3) S i2 (a,(3) //q) 

u 1 — n (1 — nj z 2 — u 
_ T ifi (a,P) T hl (a,f3) T h2 (a,(3) 
1 * a >P> u )- u + i_ n + 2-u 

and where Dij(a, /3), Sij(a, (3) and Tij(a, (3) are rational functions of a and /?. The explicit 
expressions are given in Appendix B. We note that there are simple relations between some 
of these functions. For example, for any structure function i, 

A,o(a, P) = (a + (3- 2a(3) S ifi {a, 0). (4.9) 

It is straightforward to check that there are no renormalon singularities in Bf DG (a, (3, u) 
of Eq. (4.2). As usual, single poles in bi(a, (3, u) are cancelled in B? DG (a, (3, u) by the sin(7ra) 
factor, which is associated with the gluon momentum being timelike. The double pole in 
bi(a,(3,u) at u = 1, which could have resulted in a single pole in Bf DG (a, (3, u), is in fact 
not there: according to Eqs. (4.7) and (4.8) its residue is proportional to 

(2 - k+ - K_)v / A5 i ,i(a,/3) + (k+ - K-)D iA (a,P) (4.10) 

and therefore to (a + (3)Si t i(a, (3) + A,i(ct,/3), a combination that vanishes for all i. 

Thus, Ri(a,(3) are free of infrared renormalons. Nevertheless, the series for Ri(a,(3) 
are divergent and they are Borel summable only for large enough values of (3, owing to the 
convergence constraint on the Borel integral in Eq. (4.1) for u — > 00. The consequences have 
been investigated in detail in the context of the radiative decay [17], see Sec. 2.3 there. In 
moment space, the convergence constraint is replaced by infrared renormalons [25] through 
the integration over (3 near (3 = 0, see Eq. (6.9) below. 



5. Expanding the Borel function 

In order to obtain the perturbative coeffi< 
expands the Borel function Bf DG (a, (3, u) in powers of u, see Eq. (A. 3). The expansion of 



In order to obtain the perturbative coefficients c^' \a, (3) in the second line of Eq. (4.1) one 
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2 iq^[l, 1], [2 — u],x^ is known, see e.g. Ref. [45]: 



2 Fi([l,l],[2-«],x) = -\n.{\-x)+u 



ln 2 (l-x)+Li 2 (x) 



(5.1) 



+ ir -S lt2 (x) - ln(l - x)Li 2 (x) + Li 3 (x) - - ln 3 (l - x) 
+u 3 - S 2 , 2 (x) - ln(l - x)Li 3 {x) + ln(l - x)S lj2 (x) 



+1 i n 2(i _ x )Li 2 (x) + ln 4 (l -x) + S 1>3 (x) + Li 4 (x) 



+ 



where Nielsen integrals are defined by 

\a+b-l rl 



(5.2) 



Expanding S| DG (a, /3, u) in Eq. (4.2) in powers of u, and using Eq. (5.1) to expand the 
hypergeometric functions, we obtain the coefficients, expressed in terms of the functions of 
Appendix B. At 0(a s ), using the definition of k± in Eq. (4.4) and the relation of Eq. (4.9), 
we get: 

cf\a,f5) =-A-wj -2(a + /?-a/3)ln(£) S ifi (a,p)+ (a - (3) T ifi (a, (3), J . (5.3) 



At O((3 a 2 s ) we get 



v \ ln 2 (l - k+) + Li 2 (K + ) \ ln 2 (l - /e_) + Li 2 (K_) 



K + 



K. 



■+ ! D,,, (n. D, (1 (n. J) + D,, (n. i) + ^A, 2 («,/3)^ ' ^ ' ' '" ( 1 ^ ''' " ' 



— K- 



±ln 2 (l - k+) +Li 2 ( K+ ) i ±ln 2 (l - «_) +L1 2 (k_) 



+ 



+ VA^, 1 (a,/3)-^o(a,/3) + ^ 1 (a,/3) + ^, 2 (a,/3)) (Ml^ + ML_^zT) 
+ V^T iA (a,/3) + ±T ii2 (a,/3)) 1 + Q - ln(a/3)) cfW)- 



(5.4) 



(1 2) 

Explicit expressions for c\ ' (a, (3) that are obtained upon substituting the functions in 
Appendix B and collecting terms, are listed in Appendix C. In this way it is straightforward 
to derive higher-order terms in the single-dressed-gluon approximation. 
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6. The Sudakov limit 



Having derived an explicit result for the Borel function it is straightforward to extract the 
singular terms in the (3 — ► limit, the Sudakov limit. The leading terms in this limit have 
already been computed in Ref. [25]. They have been extracted there — see Sec. 3 and 
appendix B — from an integral representation of the real-emission result for the triple- 
differential distribution, which was computed directly in the Borel formulation; here we 
re-derive it from the Borel representation (4.1) of the structure functions computed in the 
previous section based on the dispersive method. 

The Borel function is given in Eq. (4.2) where bi(a,(3,u) are explicitly written in 
Eq. (4.7). We wish to expand these results at small (3 keeping the other lightcone variable a 
as well as the Borel parameter u fixed. In this limit K- = (3(1 — a) /a + O(0 2 ) while 
k + = 1 — (3/a + O(0 2 ). This means that in Eq. (4.3) 2F1 (\\, 1], [2 — u], K-^j can be readily 

expanded at small (3 while 2-fi ^[1, 1], [2 — u], K+j cannot. To extract the leading terms at 
(3 — ► from the latter we first use the general identity: 



2 F 1 ([1,1],[2-u],k + ) = (1-u) -- 2 Fi([1,1],[1+«],1-k+) + (1- K+ 

V / U \ / Sill 7TU 

(6.1) 

The new hypergeometric function in Eq. (6.1), 2-^1 (jl,l], [1 + u],l — k + ), is of course 
expandable at (3 — ► 0, while the non-analytic contributions are explicitly given by the 
second term in Eq. (6.1). 

With this replacement and using Eqs. (4.2), (4.7) and (4.8) and the explicit expressions 
in Appendix B, we obtain the expected singularity structure [25] (see Eq. (3.17) in Ref. [16]) 
for small lightcone component (3: 



(a/3)- M ^ LO (q) e§" 
~f3 u~ 



B { \a.,p,u)\p^ = — x 



, B\ u 1 sin7ru / 1 1 

(i-u) - -o h + 



a J 2 ttu \1 — u l — u/2 



(6.2) 

x (l + 0((3/a)), 



where the last factor in Eq. (6.2) serves as a reminder that integrable terms that are 
suppressed by powers of (3 are excluded here. 

As in the full Ri(a,(3), there are additional 0(1/ do) contributions to Bi(a, (3, u)\p_> 
starting at 0(u l ), corresponding to 0(a 2 s ), the NNLO. These go beyond the large-/?o limit, 
and therefore beyond the calculation performed in the present paper. In contrast with the 
full Bi(a, (3,u), these (3 — > singular terms are known in full to the NNLO [15,16,35] - 
see e.g. Eq. (3.41) in [16] — and they play an important role in Sudakov resummation. 

The perturbative expansion in Eq. (4.1) corresponding to Eq. (6.2) contains log- 
enhanced terms of the form ~ ln k (j3) / (3 , with k < n. Going beyond the large-/3o 
limit one finds higher powers of the logarithms owing to multiple gluon emission: at the 
n-th order one obtains ln fc ((3)/ (3 where k goes up to 2n — 1, see Eq. (6.8) and Eq. (6.9) 
below. 



- 11 - 



Eq. (6.2) represents the small-/? limit of all structure functions except for i = 4 where 
the LO vanishes, so the bremsstrahlung contribution is entirely integrable, ~ ln fc (/3), 
with k < n. In the latter case we find the following leading small-/? behavior: 



2(a/3)-"e3 u 



a 



u 



1(1 - «) ( ^ 



Sill 7Tti 



ITU 



(1 



2 + a 2 + a 
+ 



u 



l-u/2 



(6.3) 



The decay width being infrared and collinear safe, the infrared singularities in Eq. (6.2) 
become integrable near the Sudakov limit once virtual corrections are included. It is there- 
fore convenient to define an integration prescription absorbing singularities from virtual 
corrections into the real emission part. Defining r = (3 /a, we use the plus prescription as 
follows: 



drF(r) 



-l+u 



F(0) 



u 



(l - r -») + jH dr(F{r) - F(0)) (6.4) 



where F(r) is a smooth test function. Upon expansion in u, this definition reproduces the 
(...)* distribution adopted in [31,32] and the plus distribution of [33]. Equivalently, we can 
use: 



1 



„1+M 



5(r) 



1 



a5{(3) 



u 



(6.5) 



Having defined the integration prescription, the real-emission coefficients c[ n \a,(3 = ra), 
at each order n in the expansion (4.1), are divided into two parts: the singular part of 
(6.2) is put under a plus prescription, and the remaining, regular part (which requires no 
prescription) is left unmodified: 



» 



+ c 



(n), 



(6.6) 



(1 2) 

Considering in particular the first two orders c\ ' (a, (3 = ra) of Eqs. (5.3) and (5.4), which 
are written explicitly in Appendix C, the part that is put under the plus prescription is: 

Vf°(a) [ ln(r) 71" 



„(l) sing. 



(a, [3 = ra) 



cf )si ^(a,(3 = ra) 



a. 



4 r 



V^(a) 



a 



31n 2 (r) 



2 r 
■ 2 85 



+ (21n(a) + - 



13 \ ln(r) 



(6.7) 



These expressions are consistent with Eq. (3.41) in Ref. [16], where at the NNLO additional 
terms, with different color factors, are included. 

The subtraction corresponding to this integration prescription will be applied when 
regularizing the virtual corrections, see Eq. (6.12) and Sec. 7 below. It should be empha- 
sized that defining plus distributions with respect to r is just a matter of convention, and 
depending of the kinematic variables used other choices may turn out more convenient. 
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In Sec. 9 we shortly discuss two alternatives, demonstrating the way in which different 
infrared-subtraction procedures shuffle terms between the real and virtual contributions. 

As discussed in Refs. [16,25,35], Eq. (6.2) contains infrared and collinear singular 
contributions that are associated with two subprocesses, which decouple in the Sudakov 
limit, the quark distribution (in an on-shell heavy quark) with momentum Oim^fi) and 
the jet with virtuality 0(mby/a/3). The two terms in the square brackets in Eq. (6.2) 
correspond to these two Sudakov anomalous dimensions, respectively, which were computed 
here, again, in the large-/3o limit. 

To perform Sudakov resummation, accounting for multiple soft and collinear radiation, 
we follow [16] and define moments of the structure functions with respect to r = (3 /a, in 
accordance with the plus prescription (6.4): 



i(a,N) = £dp(l 



N-l 



Wi(a,/?) 



Hi{a) x Sud(m fe a, A) + ARf\a), 



(6.8) 



where Hi(a) = Vi(a)\ laige ^ H . We can deduce the structure of the Sudakov exponent 

from Eq. (6.2): 



Sud(m5Ct, N) = exp < — — 




A 2 



B s (u) 



-Bj(u) 



T(N)T(-2u) 
F(N - 2u) 



1 

+ 2u~ 
1 



(6.9) 



( T(N)T(-u) 

V r(A -u) u 



which we wrote as a Borel integral (in the DGE form) with 



Bj(u) 



es l 



1 sin7rn 



large fa - 2 nu 

5„ 



1 



+ 



l-u l-u/2 



Bs(u)\ laigef3o =es u (l-u). 
In the exponent in Eq. (6.9) we added, under the Borel integral 



(6.10) 
(6.11) 



Cf r 

A) Jo 



A 2 

duT{u) ( — 2 



m 7 



the singularities that are required for writing the r 
as a plus distribution (Eq. (6.4)): 

-2u 



B[V ](a,u) 



a 



I sing. 



2u 2 



B s (u)-2Bj(u) 



non-integrable terms of Eq. (6.2) 

(6.12) 



a-^es 1 



(l-u)- 



sm iru 



7TU 



2v? 



+ 



1 



u 



1 - u/2 



25\ 1 



ln 2 (a)-|ln(a)-ivr 2 + ^) +Q{„) 



making the moments in Eqs. (6.8) and (6.9) above finite. These terms will be subtracted 
from the virtual corrections in Eq. (7.18) below. As shown explicitly in Sec. 7 this subtrac- 
tion exactly cancels the infrared singularities of the virtual corrections. 
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7. Virtual corrections 



The virtual contribution to the structure functions 1 can be decomposed as in Eq. (2.9): 

V^(p, q) = -Vi(a) 5 ^ + V 2 (a)v^v u + iV 3 (a)e^P a v p q a + V 4 (a)q p q u + F 5 (a)(/f + v u q»). 

(7.1) 

For each structure function one has an expansion in the coupling, see Eq. (7.4) below; 
according to (2.13), at leading order (Born level) we have: 

V£q( P , q) = -ag^ + Av»v v + i2e^ pc r v p q a - 2(v^q v + v v qf l ). (7.2) 

Let us now compute the virtual corrections in the large-/?o limit. To this end we modify 
the gluon propagator according to 

—k 2 ~^ {-k 2 Y+ u- { ' 

With this modification 2 , the momentum integration should yield directly the Borel function 
B[Vi\(a,u) in 

Vr°(a) - V^°{a) + f°duT(u) (^X B[^](a,u) (7.4) 

Po JO \ m b/ 

However, in contrast with the real-emission result of Eq. (4.1) that is regular for u — > 0, 
the Borel integral of the virtual diagrams is obstructed by a double pole of B[Vi](a,u) at 
u — ► 0, which corresponds to the usual double-logarithmic infrared singularity. Therefore, 
after computing the momentum integral we shall perform infrared subtraction using the 
singular terms (6.12). This would finally yield a meaningful Borel representation for the 
virtual contribution, Eq. (7.21) below. 

We define z = q 2 /m%, and in the following, since (3 = 0, we have z = I — a. The result 
of the virtual diagrams, where the gluon propagator is modified according to (7.3), takes 
the form 3 : 

B[Vn(P,q) = B[V ](a,u) x V&(p,q) 

-B[Fi](a,«)<r 
+ B[V 2 }(a,u)v»v» 

+ iB[V 3 ](a,u)e^v p q a 

+ J B[y 4 ](a,u)g /1 &, 

+ B[V 5 ](a,u) (/f + v v q?) 



1 Ova Lorentz decomposition is similar to that of Appendix B in Ref. [33]; it differs from that of Ref. [31]. 
2 See [40] or Sec. 2.2 in [44]. 

3 This result agrees with Eqs. of (B.ll) and (B.12) in Ref. [33] (with m c = 0). There a gluon mass is 
introduced instead of a Borel parameter. 
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where 



B[V 2 ](a,u) 
B[V 3 ](a,u) 



B[V }(a,u) 



B[Vi](a,u) 



ei" ^D(u) + ^C(z,u) + (l-z)(l x (z,u)+I y (z,u)-h(z,uj) ; 

-eh(l-z)(K(z,u)-I x (z,u))- 

4ei u zl xy (z,u); 

-2 el" (k(z,u)-I x {z,u))- 



B[V 5 }(a,u) 



B[V 4 ](a,u) 




(7.6) 



In Eq. (7.6) the term D{u) 



D{u) = [ dx(l - x) l+u x~ 




Jo 




(7.7) 



is the result of the 6-quark self-energy diagram' 



, while the ■u-quark self-energy diagram 



vanishes in the Borel regularization, as the momentum integral has no scale. All the 
remaining terms in (7.6) arise from the vertex correction diagram. Let us recall that the 
Borel parameter regularizes both ultraviolet and infrared logarithmic singularities, just as 
in dimensional regularization. Thus, u — > singularities in individual diagrams arise from 
both the ultraviolet and the infrared and no distinction is made between them. However, 
in the present context we know in advance that the ultraviolet divergencies cancel out in 
the sum of all diagrams — the current is conserved — and therefore the remaining u — > 
singularities in the sum of diagrams are immediately identified as infrared ones. We will 
address these singularities below. 

Let us now briefly describe the calculation of the vertex diagram and define the integrals 
entering Eq. (7.6). Upon combining the propagators using Feynman parametrization, where 
the 6-quark propagator is associated with the Feynman parameter x and the u-quark 
propagator with y (so the gluon with 1 — x — y), one identifies the scale 



where z = q 2 /m% as above. Performing next the loop-momentum integral in four dimen- 
sions one obtains integrals of the following form over the Feynman parameters: 



where a, b and c are non-negative integers. This integral is computed as follows. One 
first changes variables into from y into w = 1 — x — y and then from x into t where 

4 Our result for the self-energy diagram agrees with that of Ref. [23]. 




(7.8) 




(7.9) 
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x = (1 — w)(l — t). The integration over both w and t extends over the interval [0, 1]. 
In these variables the integrand factorizes and, assuming that the parameters a, b and c 
are such that both integrals exist (which is always the case for the required integrals) the 
result is: 

T(a-c-u + l)r(6 + l)r(l + u)r(2 + a - 2c - 2u + b) 
aAc(Z ' U) ~ r(2 + a - c - u + 6)r(3 + a + b - 2c - u) X 

2 F 1 ([c + u, b + 1], [2 + a - c - u + b], z). (7.10) 

At the next step, known hypergeometric identities (see e.g. [45]) are used to bring the result 
into one that is convenient to expand at small u. In all cases it is possible to write the 
result in terms of a single hypergeometric function 2-^1 ([1, 1 + it], [2 — it], z). This function 
has the following expansion (type E in Ref. [45]): 

2 Fi([l, 1 + u], [2 - u], z) ~ I - ln(l - z) - u ( - ln 2 (l - z) - Li 2 (z)) 

+ u 2 ( - 2S lt2 (z) - 21n(l - z)U 2 {z) + Li 3 (z) - | l n 3 (l - z)) 

- (2,S2,2(^) + 21n(l - z)Lis(z) - 41n(l - z)5i, 2 (^) 

- 2 ln 2 (l - z)Li 2 (z) - ^ ln 4 (l - z) - 4Si, 3 (*) - Li 4 (z)) + 0{u A ) | . 

(7.11) 

Following Ref. [33] (see Eq. (B.3) there) we separate the numerator of the vertex 
diagram into N± u , which is composed of terms having powers of the loop momentum in 
the numerator (cf. Eq. (B.6) in Ref. [33]) and other terms, N% v (cf. Eq. (B.9) in Ref. [33]). 
In Eq. (7.6) above, the N^ u gives rise to C(z,u) in B[Vo](a,u): 

C(z,u) = ^K(z,u), (7.12) 

where 

"1 /-l-x ^2ui 



r(i-2u)r(i + «) 



r(3 - u) 



l-^-^-2F 1 ([l,l + u],[2-u],z) 
1 — u 



(7.13) 



i + „Q + L_iM(i-z)) + fV). 



while the iVg 4 " is the source of all the other terms, where the following additional integrals 
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show up: 
h(z,u) 



1 , f X ~ x , m 2u (l-x-yY 
ax / ay- 
Jo 



(M 2 ) 1+u 



1 + 



r(i + u)r(i -2u) 

2u 2 T{l-u){\- z) 



1 ] 1 ln(l-z) 



2uz 
1 - it 



2 F 1 ([l,l + u],[2-u],z) 



l-z] 2u 2 



u 



+ ^ + ln 2 (l - z) + Li 2 (z) - | ( ln(l - z)vr 2 
- 3C 3 + 6Si, 2 (*) +61n(l - z)Li 2 (z) - 3Li 3 (z) + 21n 3 (l - z)) + 0(u 2 )|, 



= [ dx [ 
Jo Jo 



(7.14) 



I x (z,u) 



(M 2 ) 



r(i - 2u)r(i +«) 



2 Fi([l,l + u],[2-u],z) 



(i-u)r(2-u) 

= - ln(1 ~ Z) + i (ln 2 (l - z) + Li 2 (z) - ln(l - z)) u + 0(u 2 ), 



(7.15) 



Iy{z,u) 



fir/"-"''' 

JO jo 



(M 2 ) 



r(i + u)r(i - 2u) 

' uT{2-u){\- z) 



1 + ! T^r 2Fl([1 ' 1 + nU2 " u] ' z) 



j^—l ---l+(l + -)ln(l-z) 
1 — z I u \ z / 



+ 



1 + - )(- ln 2 (l - z) + ln(l - z) - Li 2 (z)) - ^vr 2 - 1 

2^ / O 



and 

Ixy(z 



Jo Jo 



io Jo (M 2 ) 

r(i + u)r(i - 2u) 



zr(3 - u) 



l+u 
l-(l + Z 



1 - u 



2 Fi([l,l + u],[2-u],z) 



(7.16) 



u +0(n 2 ) L 



2z 2 ' 



ln(l - z) - z + 
u + 0(u 2 )\ 



ln 2 (l-z)+ z-- ln(l-z) 



(7.17) 



- -z + U 2 (z) 



Eq. (7.5) is written such that all the u — > infrared singularities are in the first term. 
This term is proportional to the LO result, as must be the case. As mentioned above, in the 
absence of such singularities one would interpret the expansion of the virtual corrections, 
starting at 0(a s ), according to Eq. (7.4) above. Obviously, since there is a double pole 
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at u = 0, the u-integral (7.4) is ill-defined. One should first perform subtraction of the 
singularities. 

Using the explicit results for the integrals given above, it is straightforward to check 
that the u — ► singularities of B[Vo](a,u) in Eq. (7.6) do indeed coincide with Eq. (6.12), 
that was determined in Sec. 6 by defining the plus distribution for the real-emission terms. 
After subtracting from B[Vo](a,u) in Eq. (7.6) the terms in Eq. (6.12) one has 



B[V^](p,q)\ 



B[V }(a,u)- B[V ](a,u)\ sing . 

- B\yi](a,u) g^ v 
+ B[V 2 ](a,u) v^v v 
+ iB[V 3 ](a,u) e^v p q a 
+ BlV^a.ujq^ 
+ B[V 5 ](a,u) + 



x V£(p,q) 



(7.18) 



Next, let us split the regularized terms proportional to V^q using Eq. (7.2) and absorb 
them into the five different structure functions; we define: 



B[H(M)| reg .=-5[fi](a,u)| reg . sT 
+ B[V 2 ](a,u)\ IC& v»v v 
+ iB[V 3 ](a,u)\ ieg . e^v p q a 
+ B[V 4 ]{a,u)\ reg q^q u 
+ B[V 5 ](a,u)\ (v»r + v"e). 



(7.19) 



Finally, using Eq. (7.6) with the explicit results for the integrals we get: 



| J (2q + «-2) 



a 



-2u 



2v? 



1-u 



u-l 

sin ttu 



G(u) ¥(u, a) + ( - + 1 + \u - \u 2 ) G(u) 



1 1 

+ 



TTU \ 1 — U 1 — u/2 



B[V 2 ](a, u)| reg . = - S[Ki](a, u)| reg . + 1 }J > G(«) F(u, a) 



B[y 3 ](«,«)l rcg . = - B[Fi](a,« 



a 



ircg. 



4-7/ r > 

5M(«,n)| reg .~e^G(,) 



(2a - u - 1) 
u-l 



F(u,a) + 1 



2 2?/p3 u 
B[F 5 ](a,u)| reg . = -- B[Fi](a,u)| reg . + — -G(«) 



2(u + l) -a{u + 3) 
u-l 



¥(u,a)-l , 
(7.20) 



where G(u) = ~r(l - 2tt)r(u)/r(3 - it) andF(u,a) = 2 -Fi([l, l + [2-u], 1 - a) 
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The infrared-subtracted virtual terms can be expanded order by order in u to get 
the perturbative corrections in the large— po limit. Let us write, in analogy with the real- 
emission result (4.1), a Borel representation of Vi(a) in (2.12): 



Vr G ( a ) = V^(a) + ±f duT( 

Po Jo 

= V l h °{a) + C F 



( A 2 



2 \ u 



B[VA(a,u)\ 



reg. 



(7.21) 



(i) a s {m b ) ( 2 ) 



TT 



a s (m b ) 



TT 



+ 



fl 2) 

The coefficients ' (a) for the five structure functions, i = 1 to 5 are listed in Appendix D. 



8. The triple— differential width at NLLO in the large— f3 limit 

In the previous sections we computed separately real and virtual corrections to the five 
different structure functions in the decomposition of the hadronic tensor. For massless lep- 
tons, only three of those enter into the expression for the spectrum (2.6) through (2.11). Let 
us now combine the result into an expression for the triple-differential width. We present 
explicit expressions up to NNLO which are valid in the on-shell quark mass scheme. Using 
the results of the previous sections the generalization to higher orders is straightforward. 

Let us write the perturbative expansion of the triple-differential width in the large-po 
limit as follows: 

2 . (8.1) 



+ C F 



ot s (m b ) (a s (m h )\ 



IT \ IT 



where a>o(a, x) is defined in (2.14) and the NLO and NNLO coefficients K n (a, (3, x) for n = 1 
and 2, respectively, will be detailed below. At each order real and virtual contributions to 
each structure function add up according to (2.12). The coefficients of the triple differential 
width can therefore be written as follows: 

K n (a,p,x) =Lo n {a,x)5{(3) + {K^{a,(3,x)} + + K^ (a, (3, x), (8.2) 

where, as usual (3 = ar and the plus prescription is as defined with respect to r according 
to Eq. (6.4). 

To obtain the virtual coefficients u>i(a,x) at each order n one substitutes the regu- 
larized virtual coefficients of the structure functions (7.21), which are given explicitly in 
Appendix D for n = 1, 2, into Eq. (2.11) and uses the result in (2.6). The virtual coefficient 
at NLO (n = 1) is 

( 5 7T 2 \ 

wi(a,x) = - u (a,x)l Li 2 (l - a) + - + — J + 6 (a - 1 + x) (2 a - 5 + 2 x) In (a) (8.3) 
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and the NNLO result (n = 2) is: 



0^2(0, x) = 6 (a — 1 + x 
- u (a,x) 



25 

Li 2 (l - a) + ln 2 (a) - — ln(a) 



6 



19. 



Li 3 (l - a) + 2Li 3 (a) + — Li 2 (l - a) 

o 

+ (ln(l - a) - 1) ln 2 (a) + _ - - tt 2 J ln(a) - ( 3 + — + — 



•4) 



Similarly, by substituting the real-emission results for the structure functions into 
(2.11), one obtains the corresponding real-emission coefficients for the triple differential 
width. The singular part, K^ ng -(a, (3,x), which enters (8.2) under a plus prescription (6.4), 
is obtained by 

K^ ne '(a, ar,x) =6(1 - a)(l - ar) smg '(a, a r) 

— 3 ^x 2 — x (2 — a — a r) + (1 — a) (1 — a r)^j c 2 "^ sms '(a, a r) ^ ^ 

+ 6(l-a)(l-ar) (x - 1 + i(a + ar)J cj sing '(a, a r) . 

Since at any order n the coefficients c^ sins '(a,ar) are proportional to the corresponding 
LO result V^°(a) for i = 1 to 3, one obtains the singular part i^ ins '(a, /3, cc) as the r — ► 
singular terms corresponding to the expansion of (6.2) in powers of u, which depend only 
on r and a, times the following prefactor, 



Q(a, r, x) = d>o(a, x) + 6^2a 2 — 7a — 4x + 5 + 2xaj ar — 6a 2 (l — a) r 2 , 



(8.6) 



that depends also on the lepton energy fraction x. In particular, using Eq. (6.7), the NLO 
result is 



i^ mg '(a, ar, x) 



Sl(a, r, x) 



a 



and the NNLO one is: 



K^'^a, ar,x) = 



fl(a, r, x) 



a 



31n 2 (r) . 



ln(r) 7 1 
r 4 r 



13 \ ln(r) (7 



(8.7) 



vr 2 85\ 1 



+ -Ma) + 
r V 2 V 7 6 24 / r 



As expected, the O(r ) term in Q,(a, r, x) coincides with the Born-level result uo(a,x). 
However, owing to the contraction with the leptonic tensor in (8.5), Cl(a, r, x) also con- 
tains some ©(r 1 ) and 0(r 2 ) terms that generate integrable O(r ) and 0(r l ) terms in 
K^ ng (a,ar,x). These terms can be freely taken out of the {...} + brackets in (8.2), as they 
do not vary by applying the plus prescription (6.4). 
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Finally, the regular terms K™ g '(a, (3,x) are obtained by 
K™ g '{a, a r, x) =6(1 - a)(l - ar) c["'(a,ar) — c^ miie '(a,ar) 



3 ( x 2 — x (2 — a — a r) + (1 — a)(l — ar) ) df 1 ' (a,ar) — c 2 "^ smg '(a, a r 



+ 6 (1 — a)(l — ar) — 1 + -(a + ar)J c^\a,ar) — c^ s ™ e (a, a r) 

where, for n = 1 and 2, the explicit expressions for c^\a, (3) and its singular part are 
given in Appendix C and in Eq. (6.7), respectively. Using these expressions we find that 
the regular term at NLO is given by: 



(8.9) 



K^ (a,ar,x) =J^lQiK r, x) 



a(l -r) 4 

3(l-a)(l-ar) 
2(l-r) 



3(1 — x — a) (1 — x — cir) . ,„ 

> 2 (a, r, x) - -± {1 _ r) 3 a r, x) (8.10) 



3(l-a)(l-ar)(2x-2 + a + ar) , 

'-(4a 2 r- 10ar + 7r- 10 a + 9), 



2a(l -r) 



and at NNLO by 



K T 2 ce '(a, a r, x) = — 2KJ CE '(a,ar, x) ln(a) — 3 



1 ln(r)P2(a, ^ 9Qi(a, r, x) 2 
+ ~ t: ; — \T r; In (r) + 



(1 + r — a r) ln(l + r — a r) 
ttr(l-a)(l-ar)(l-r)3 r ' X) 

6Qi(a, r, x) f2(a, r, x) 



+ 
+ 



2 (l-ar)a(l-r) 4 a(l-r) 4 

r (1 — a) 
1 + r — a r 

1 P3(a, r, i) 1 (1 — x — a r) (1 — x — a) 



ln(l + r — ar) ln(r) + Li 2 
+ 



a(l -r) 4 
1 — a r 



a r 



4 (1-r) ' 2 (l-r) 3 a 
1 (2 x - 2 + a + a r) (1 - a r) (1 - a) 
4 



— Ll2 . 

1 + r — a r 
F^(a, r, x) 

7r 2 Q(a, r, x) 



(1 — r) a 



Ps(a, r, x 



6 a r 



.11 



where f2(a, r, x) is given in Eq. (8.6) and the polynomials Qj(a, r, x) for j = 1 to 3 and 
Pj(a, r, x) for j = 1 to 5 are listed in Appendix E. 



9. Changing variables: alternative subtraction procedures 

In the previous sections we have presented the results for the triple-differential b — ► ulP 
width to all orders in the large-/?o limit. We have chosen to describe the hadronic tensor 
in terms of the lightcone variables a and (3 and defined plus distributions with respect to 
r = (5/a. Let us now shortly describe how the results can be used with other kinematic 
variables. This is often useful for deriving analytic expressions for partially-integrated 
spectra, as done for example in Ref. [31] at the NLO level. 

While the result for the virtual diagrams in a given regularization is unique — in 
the Borel regularization it is given by Eqs. (7.6) and (7.5) — the infrared subtraction 
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that renders their perturbative expansion finite, namely Eqs. (7.18) and (7.19), crucially 
depends on the corresponding real-emission terms that are put under the plus prescription, 
Eq. (6.2), and the variable (r) with respect to which the plus prescription is defined, see 
Eqs. (6.4) and (6.12). 

Let us first demonstrate how to use the results of the previous sections in the same 
kinematic variables a and (3, but with a different infrared-subtraction convention. Consider 
defining the plus distributions with respect to (3; instead of Eq. (6.4) we now write 



d/3F(J3) 



o 



f3^ 



F(0) 



u 



(l-/3 - u ) + j\(3[F((3)-F(0))^, (9.1) 



that corresponds to the replacement 

1 



6<J3) 



u 



(9.2) 



Taking the corresponding moments and applying (9.1) we get (cf. Eq. (6.8)): 



Wf\a,u) = / ^(l-Z^W; 
Jo 



(a, (3) 



H { /\a) x Sud ((S) (m b ,a,v) + Ai?f } (a), 



(9.3) 



where the superscript ((3) is used to distinguish the current definition from our default one, 
and 



Sud^(mb, a, u) = exp < ^f- [ 

I Po Jo u 



du 



T{u) 



mi 



B s (u) 



T(u)T(-2u) 1 
T(u - 2u) 2u 



a- u Bj(u) 



rHr(-n) | i 
r(i/ — u) u 



(9.4) 



where the large-/?o anomalous dimensions Bs{u) and Bj(u) are given in (6.10). Note that 
the explicit a dependence in (9.4) is different from (6.9) owing to the different meaning of 
the moment variable v compared to N. The subtraction term B[Vo](a,u)\ siiig ^, replacing 
(6.12), is therefore: 



B[V ](a,u)\ singm = ^ 



B s (u)-2a~ u Bj(u) 



(9.5) 



es" 
2^ 

1 

~2^2 



(1 - u) - a- u 



sm iru 



+ ln(a) - 



TTU 

25\ 1 



1 1 

+ 



12 J u 



1 — it l-u/2 

1 . o, , 29 

2 "" v " y 12 



, 2/ x ^, / n 1 2 245 

ln 2 (a)-— ln(a)--7r 2 + — 



Finally, using (9.5) in Eq. (7.18) we get the corresponding infrared-subtracted version of 
the virtual terms that replaces (7.19) in this alternative convention. The final results for 
the virtual terms, equivalent to (7.20) immediately follow. 



- 22 - 



In a similar way one can consider the subtraction using other kinematic variables. A 
natural choice, which was used in [33] as well as in the original derivation of the singular 
terms in [25], is based on the invariant masses of the hadronic and the leptonic systems, 
Eq. (2.4). Let us define: 



1-Z = % = a(3; 
wP(z,0=W i (a,/3) 



z = \ = {\-a){\-ff) 
m b 



d(a, (3) 



1 



(9.6) 



a — (3 



Here, infrared singularities are associated with the small jet mass limit — > 0, 



Brw>«)i € ^i 
(i-«) 



V^°(a = 1-z) ei l 
(1 - C) 1+M u 
1 — £ \ ~ u 1 siniru 



(1 - zf 



1 1 

+ 



2 iru \ 1 — u 1 — u/2 



(i + o(i-0) 



(9.7) 



so plus distributions are defined with respect to £ (i.e. subtracting 5(1 — £) terms) and the 
corresponding moments are: 



t(0, 



z,n)= I <Z££ ; 
o 

r(0 



2,6 



#P(z) x Sud^(m b ,z,n) + Ai?P(z), 



(9.8) 



with 



Sud«W,n) =exp(^ r-r(u) 



r(n)r(-u) | l 

r(n — u) u 



(9.9) 



Therefore, in these variables the subtraction term takes the form 



B[V }(z,u)\ singm ee — 2 



(l-zf u B s (u) - 2Bj(u) 



es" 
2^2 



(1 - z) 2u (1 - u) 



sin 7rti 



1 1 

+ 



7TU \ 1 — U 1 — u/2 



(910) 



_L + ln(1 -,)-§)!-(- >»*(! - ,) - 1 Mi - ,) - y + f | + ot ,„. 



As before, the corresponding infrared-subtracted virtual terms can be obtained using (9.10) 
in Eqs. (7.18) with a — > 1 — z. 
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10. Conclusions 



We have computed the perturbative expansion of the triple differential width in b — ► X U W, 
to all orders in the large-/3o limit. This is an important step in determining the differential 
spectrum beyond the NLO. 

Several independent partial calculations have been done in the past that provided 
useful checks. We find complete agreement with the following: 

• The NLO calculation of the five structure functions (or the fully differential width) 
in Ref. [31]. 

• The NNLO result of Ref. [21], Eq. (1.1) above, where we could check the (3q piece 
upon performing phase-space integration according to (2.7). 

• The NNLO single-differential distribution with respect to p+, computed in the large- 
Po limit in Ref. [32], which we checked by defining the subtraction procedure with 
respect to (3 = l m b (see Sec. 9) and then integrating over x and a. 

• The singularity structure of the real-emission terms as a function of the Borel variable 
in [25], and the corresponding Sudakov exponent [16]. 

• The results of Ref. [33], which have been used here for the real-emission diagrams 
and computed by a different method for the virtual ones. 

As explained in the introduction, it is expected that the O{(3oa 2 s ) contribution com- 
puted here constitutes the bulk of the 0(a 2 s ) correction. It therefore has an immediate 
application in improving the calculation of partial branching fractions used in the determi- 
nation of | Vub | from inclusive measurements in the B-factories with a variety of kinematic 
cuts. 

Although higher-order 0(/?Q _1 ag) corrections, n > 3, may also be significant, we do 
not expect that direct use of the single-dressed-gluon result by itself, namely (4.1) and 
(7.21) would yield a viable description of the triple differential spectrum. Owing to the 
u — > oo convergence constraint, the Borel integral of the real-emission corrections (4.1) 
does not exist 5 for small p~j~, namely in the Sudakov region. Better treatment of this region 
is achieved using moment space [16,25], where the convergence constraint is replaced by 
infrared renormalons. Moreover, in the Sudakov region the effect of multiple soft and 
collinear radiation is very important, and can be taken into account by exponentiation 
(6.8), as done in Ref. [16]. The running-coupling corrections computed here can be used 
to improve the calculation of partial branching fractions in the DGE approach of Ref. [16] 
by incorporating the residual O{(3oa 2 s ) correction that is not part of the Sudakov factor 
into the matching coefficient. 

Higher-order running-coupling corrections are also useful for understanding the inter- 
play between perturbative and non-perturbative corrections, and for estimating the latter. 
Our final results for the hadronic tensor, written as analytic functions in the Borel plane, 

5 The implications have been studied in detail in Ref. [17] in the context of b — > X s ^\ see Sec. 2.3 there. 
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can be used to determine infrared renormalon ambiguities, which are in turn indicative of 
the form and potential magnitude of power corrections. In the virtual part of (7.21) with 
(7.20) one identifies infrared renormalon singularities at integer and half integer values of u. 
In contrast, the differential real-emission contribution, of (4.1) with (4.2), does not present 
renormalon singularities; as mentioned above these do show up in moment space owing 
to the integration over p^ near the singular pj — > limit. At the level of the Sudakov 
exponent this has already been observed in Ref. [25] and been put to use in Ref. [16]. The 
present results facilitate the analysis of power corrections over the entire phase space. 

Let us end with a brief comment on the technical tools used in this paper, which 
made it possible to derive analytic expressions for the Borel transform. We exploited two 
different techniques for the calculation of Feynman diagrams with a single dressed gluon: 

• The real-emission diagrams were computed in Ref. [33] using the dispersive approach, 
where the gluon in the final state is assigned a fixed virtuality, which is then used 
in a dispersive integral with the time-like discontinuity of the coupling. Here we 
converted the result into a Borel representation and derived analytic expressions for 
the Borel function. 

• The virtual diagrams were computed here directly in terms of the Borel variable, 
by modifying the gluon propagator according to Eq. (7.3), and then preforming the 
momentum integral. 

Having brought the results of both real and virtual diagrams with a single dressed gluon 
to a common regularization, we could directly perform an all-order infrared subtraction. 
In this novel approach the Borel variable has a double role: on the one hand it serves 
as an infrared regulator for logarithmic singularities — a double pole at u = 0, in full 
analogy with dimensional regularization — and on the other, it serves as a conjugate to 
mm 2 /A 2 , or the inverse of the coupling constant, allowing for all-order resummation of 
running-coupling corrections. 

Acknowledgements 

The work of PG is supported in part by MIUR under contract 2004021808-009. 
A. The scheme-invariant Borel transform 

In this paper, e.g. in Eqs. (4.1) and (7.21), we use the scheme-invariant [46] Borel repre- 
sentation where T(u) is defined as the inverse Laplace transform of the coupling: 



Here the Borel variable u is the Laplace conjugate to the logarithm In// 2 /A 2 , rather than 
to the inverse of the coupling constant, which is used in the standard Borel transform. In 
this way it is possible to promote the calculation performed in the large-/?o limit to include 




(A.l) 
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running-coupling terms that are associated with subleading corrections to the (3 function, 
in particular f3\. This was used in various applications, see e.g. Ref. [47]. 

In the strict large-/?o limit, one resums the terms /?Q _1 a™ to any n. In this case 
T(u) = 1, and Eq. (4.1) (or Eq. (7.21) for the virtual terms) reduces to the standard Borel 
transform, with respect to ^ms(^) = Poaf s (/j,)/7r, namely 

i?! arsc/3 >,/3) = ^f r du exp{ - u/A m (jj,)\ B™ G (a,P,u) (A.2) 

Po jo 1 } \ m bJ 

n \ We a\ a s( m b) . (2), a , a ( a s (m b )\ 2 

TT \ TT J 



with the relation: 



CO 



5f G (a,/3,n) = ^cf t) (a,/3)^ T . (A.3) 



n 

n=l 



Upon using the scheme-invariant formulation as in Eq. (4.1), it is straightforward to include 
Pi effects in the running of the coupling by introducing T(u) that corresponds to the Laplace 
transform of the two-loop coupling (or the 't Hooft-scheme coupling): 

s f3 a* u °° ft (p,) f°° , m , , fA 2 \ u dA a2 ^ _ a . 

= \ =L duT{u) \^) ; d^=- A{i+5A) ^ 
r <«> = w£r ln( ^ A2) 4- 5in ( i+ A)' < A - 4 > 

with 5 = Pi/Pq. Upon expanding Eq. (4.1), or Eq. (7.21), with T(u) of (A.4) one obtains, 
in addition to the large-/3o terms, also P\Pq~ 3 o^ terms, etc. 

B. The functions of the lightcone variables entering the Borel transform 

Here we give explicit results for the rational functions of the lightcone variables, Dij(a, P), 
Si j (a, P) and Tij(a, P) entering the Borel function of Eq. (4.8). Recall that the first index 
here, i = 1 to 5, corresponds to the structure function, while the second to the location of 
the singularity on the positive real axis in the Borel plane, u = j with j = to 2. There 
are a few relations among them, valid for any i: 

D ifi (a, P) = (a + P-2aP) S ifi (a, P) (B.l) 
D i: i(a,P) = -(a + P)S i: i(a,P) (B.2) 



D it 2{a, P) = ~ — 2 2 ; ^ a/ , ^ 5 i)2 (a, P) (B.3) 



a 3 (l- p)+ P 3 (l-a) 
a 2 + P 2 + aP(l -a-P) 

The remaining functions are: 



(-a- p + a 2 + p 2 + 3ap + 3a 2 P 2 -3a 2 P-3aP 2 ) 
1,2(a '^" 2(-l + p) (-1 + a) (-p - a + a p) X (B.4) 

(-a 2 + a 2 p-aP + ap 2 - p 2 ) 



- 26 - 



Su(a,0) =(-a 2 + 2a 2 /3 3 + 2a 3 -/3 2 -3a 3 /3 + 2/3 3 + 2a 2 /? + 2a 3 /3 2 -4a 2 /3 2 
+ 2a/3 - 3a/3 3 + 2a/3 2 )/(2 (-/? - a + a/3)) 

51,1(0!, /3) = - a(- 11 a 3 /3 + 5 a 3 + 6 a 3 /3 2 + 31 a 2 p - 15 a 2 - 22 a 2 /3 2 + 6 a 2 /? 3 

- 34 a (3 + 12 a - 11 a /3 3 + 31 a (3 2 - 15 /3 2 + 5 p 3 + 12 0)/ (B.6) 
(4 (-1 + 0) (-1 + a)(-/3 - a + a/3)) 



2/? 4 - 2a/3 4 + a 2 /3 4 - 4a z P A - 4a z p' z - 4a 3 Z + a 4 /?^ - 2a 4 p + 2 a 



51,0(0, P) 



v2 «3 



2 fl2 



v 3 fl2 , „,4 «2 



4a/3(-/3-a + a/3) 



(B.7) 



D^! ( a , /3) =(3 a 4 /3 3 + 13 a 4 /3 + 3 a 3 /3 4 - 16 a 2 /3 4 + 13 a /3 4 - 16 a 4 /3 2 + 2 a 3 /3 5 + 2 a 5 /? 3 
+ a/3 5 + a 5 /3 - 3a 5 p 2 - 3a 2 /3 5 + 4a 3 + 11 a 3 p 2 - 8a 3 /3 + 11a 2 /3 3 
+ 8 a 2 /3 2 - 4 a 2 P - 8a p 3 - 4a p 2 + 4/3 3 - 10 a 3 p 3 + 4 a 4 /3 4 - 4a 4 
- 4/? 4 )/(4 (-1 + P){-\ + a) (-/? - a + a/3)) 

(B.8) 

T 1:2 (a,p) =- (6a 4 /3 3 - 8a 4 /3 2 + a 4 /3 + a 4 + 9 a 3 /3 2 - 2a 3 /3- a 3 + 6a 3 /3 4 - 8 a 3 /? 3 
- 6 a 2 /3 2 + a 2 /3 + 9 a 2 /3 3 - 8 a 2 /? 4 + a/3 4 - 2a/3 3 + a p 2 - p 3 + /3 4 )/ 
(4(-l + /3)(-l + a)a/3) 

(B.9) 

Ti,i(a,/3) =(4a 4 /3 3 - 8 a 4 /3 2 + 5 a 4 /3 - a 4 - 12 a 3 /3 3 + 13 a 3 /3 2 - 4 a 3 p + a 3 + 4 a 3 /3 4 
-8a 2 /3 4 -6a 2 /3 2 -a 2 /3 + 13a 2 /3 3 + 5 a /3 4 - 4 a /3 3 - a/3 2 + /3 3 - /3 4 )/(2 
(-l + /3)(-l + a)a/3) 

(B.10) 



Ti, {a,P) =- 



-7aP(a + P) - Wa 3 p -10aP 3 + 2a 2 p 3 + 7/3 3 -4a 2 /3 2 + 2 a 3 p 2 + 7 a 3 



4ap 



(B.11) 



52,2(0, /3) = - 2 (9 a 2 /3 2 - 9 a 2 p + a 2 - 9 a /3 2 + 13 a p - 3 a + /3 2 - 3 /?) 
(-a 2 + a 2 p-ap + aP 2 - P 2 )/(~P - a + a/3) 



5 2 ,i (a, P) = - 2 (9 a /3 2 - 2 /? 4 + 39 a 2 /3 3 + 6 a 4 /3 3 + 39 a 3 /3 2 - 34 a 2 /3 2 + 9 a 2 /3 - 2 a 4 
+ 3a 3 - 20a 3 /3- 13a 4 /3 2 - 20a/3 3 + 9a/3 4 + 9a 4 /3- 28a 3 /? 3 - 13a 2 /? 4 
+ 6 a 3 /3 4 + 3 p 3 )/{-p -a + ap) 

(B.13) 
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S 2 ,i(a, (3) =(-35 a 4 (3 2 + 18 a 4 f' • 17n ' i • 91 a 3 (3 2 - 62 a 3 /3 3 + 6 a 3 - 52 a 3 /? 
/3 4 -35< 

+ 30 a /3 2 + 6 /3 3 )/(-/3 -a + a/3) 



+ 18 a 3 /3 4 - 35 a 2 (3 4 - 88 a 2 f3 2 + 30 a 2 /? + 91 a 2 /3 3 + 17 a p 4 - 52 a /3 3 



(B.14) 

S 2 ,o{a, P)=- (-21 a 4 /? 3 - 2 a 4 /3 - 21 a 3 /3 4 + 18 a 2 /3 4 - 2 a /3 4 + 34 a 3 /? 3 + 4 a 4 p 4 
- 14 a 3 /? 2 - 14 a 2 f3 3 + a 5 (3 3 + 4 a /3 5 + 18 a 4 /3 2 + a 3 /3 5 - 3 a 2 /3 5 - 2 a 5 
+ 4a 5 /3-3a 5 /? 2 - 2 /3 5 )/(a (-/? - a + a (3)) 

(B.15) 



D 2jl (a, /?) = - (- 12 a f3 3 - 12 a 3 + 27 a 3 /3 4 + 14 a 5 /3 3 + 14 a 3 /3 5 - 27 a 2 /3 5 - 27 a 5 /3 2 
- 4a/3 4 - 20a 2 /3 4 - 20 a 4 /3 2 + 2 p 4 + 2 a 4 - 4a 4 + 27a 4 /3 3 + 76a 3 [3 2 
+ 76a 2 /3 3 + 13a/3 5 + 13a 5 /?- 106a 3 /3 3 -4a 4 /3 4 
-28 a 2 (3 2 )/{-(3- a + a/3) 

(B.16) 

T 2j2 (a, /?) =(18 a 4 (3 3 - 24 a 4 /? 2 + 5 a 4 /? + a 4 + 25 a 3 /? 2 - 14 a 3 /? + 18 a 3 /3 4 - 24 a 3 /? 3 
+ 2 a 2 [3 2 + 25 a 2 /? 3 - 24 a 2 p 4 + 5 a /3 4 - 14 a /3 3 + /3 4 )/(a p") 

(B.17) 

T 2jl (a, /?) = - 2 (-/? + 2 a /3 - a) (-7 a 3 /? + 6 a 3 (3 2 + a 3 + 6 a 2 /3 3 - 16 a 2 /3 2 
+ 11 a 2 /? - 7 a /3 3 + 11 a /3 2 + /3 3 )/(a/3) 

T 2j0 (a, 13)=- (-7 p 4 - 20 a 4 /3 2 - 10 a 3 /? - 20 a 2 /3 4 - 10 a /3 3 - 56 a 3 /3 3 + 21 a 4 /? - 7 a 4 
- 38 a 2 /3 2 + 57a 3 (3 2 + 21 a /3 4 + 57a 2 /3 3 + 6 a 4 /? 3 + 6 a 3 /? 4 )/(a p") 

(B.19) 

-a 2 + a 2 /?-a/3 + a/3 2 -/3 2 

53,2(a, p)= t, : a ( B - 20 ) 

— p — a + a (3 

~ . -3a/3 2 + 2/3 2 + 2a/3 + 2a 2 /? 2 -3a 2 /3-/? + 2a 2 -a 

#3,1 («>/*) = ■ (B.21) 

— (3 — a + a/3 

(3a(2a(3-3a + A-3(3) 
53 ' l(a)/3) - 2(-/3-a + a/3) (R22) 

_ , ox -2ap 3 + a 3 /3 2 + a 2 /? 3 + 2a 3 + 2/3 3 + 2a 2 /?-6a 2 /3 2 -2a 3 /3 + 2a/? 2 

<33,o(a,P) = — — ^ ■ -rr 

2ap (— p — a + ap) 

(B.23) 
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(a 2 /? + a/3 2 -8a/3 + 4a + 4/?)(-/? + 2a/?-g) 
D v( a >® = 2(-/?-a + a/3) (R24) 

2a 2 /3 + a 2 -2a/3 + 2a/3 2 + /3 2 
T 3>2 (a, /J) = ^ (B.25) 

(-/3 + 2a/3-a) 2 
r s ,i(a,/3)= ^ (B.26) 

„ , 7a 2 + 7/3 2 + 4a 2 /3 2 -10a/3 2 + 2a/3-10a 2 /3 , . 

rs,o(a, /?) = ^ (B - 27) 

S 4)2 (a,/3) =- (-a 2 + a 2 /?-a/3 + a/3 2 -/3 2 )(4a 2 - 6 a - 6 /3 + 4/3 2 + 22 a /3 + 12 a 2 /? 2 

- 18 a 2 /3- 18a/3 2 + 3a/3 3 + 3a 3 /?)/((-l + /?) (-1 + a) (-/? - a + a /?)) 

(B.28) 

S 4 ,i(a, /3) = - 2(2 /3 3 + 5 a 2 /3 3 + 2 a 3 + 10 a 2 - 8 a /3 3 + 5 a 3 /3 2 - 14 a 2 /3 2 + 10 a f5 2 
+ a p 4 + a 4 /? - 8 a 3 - a + a/3) 

(B.29) 

5 4jl (a, /3) =(195 a 4 /3 3 + 78 a 4 /? + 195 a 3 p 4 - 187 a 2 p 4 + 78 a p 4 - 187 a 4 /? 2 - 51 a 3 /? 5 

- 51 a 5 /3 3 - 28 a /3 5 - 28 a 5 (3 + 60 a 5 /3 2 + 15 a 5 p 4 + 15 a 4 (3 5 + 60 a 2 /? 5 

+ 6a 3 + 3a 3 /3 6 + 3a 6 /3 3 + 3a 6 /3 + 3a/3 6 -6a 2 /? 6 - 6a 6 /3 2 + 4/3 5 + 4a 5 
+ 254 a 3 /3 2 - 84 a 3 /? + 254 a 2 /3 3 - 148 a 2 (3 2 + 30 a 2 p" - 84 a /3 3 + 30 a /3 2 

+ 6/? 3 - 326a 3 /? 3 -90a 4 p 4 - 10a 4 - 10/3 4 ) /((-/?- a + a /3) (-1 + a) 2 
(-1 + /3) 2 ) 

(B.30) 

, m 2(-2a 3 + a/3 3 -2/? 3 -7a 2 /3 + a 3 /3 + a 2 /3 2 -7a/3 2 ) 

5 , 4 ,o(a,p)= s . 3 (B.31) 

— p — a + ap 

D 41 ( a , /?) = - (3 a 7 /? + 207 a 4 (3 3 + 6 a 4 /3 + 207 a 3 p 4 - 116 a 2 p 4 + 6 a p 4 - 116 a 4 (3 2 

- 29 a 3 /? 5 - 29 a 5 /? 3 + 22 a /3 5 + 22 a 5 /? + 5 a 5 /? 2 + 5 a 5 p 4 + 5 a 4 /3 5 
+ 5 a 2 /3 5 + 3 a 3 /f + 6 a 6 p 4 + 6 a 4 /3 6 + 6 a 5 /3 5 + 3 a 7 /3 3 - 23 a 3 /3 6 

- 23 a 6 (3 3 - 20 a 6 /? - 20 a /3 6 + 33 a 2 /3 6 + 33 a 6 /3 2 + 4 a 6 + 4 /3 6 
+ 3 a /3 7 - 6 a 2 (3 7 - 10 /? 5 - 10 a 5 - 6 a 7 /? 2 + 124 a 3 p" 2 - 12 a 3 /? 
+ 124a 2 (3 3 - 36a 2 (3 2 - 12a /3 3 - 276a 3 /? 3 

- 110a 4 /? 4 + 6a 4 + 6/? 4 ) /((-/? - a + a /3) (-1 + a) 2 (-1 + /3) 2 ) 

(B.32) 
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T 4)2 (a, 0) =(3 a 4 - a 4 + 15 a 3 2 - 20 a 3 /3 + 7 a 3 - 30 a 2 2 + 23 a 2 - 10 a 2 + 15 a 2 /3 3 
+ 3a/3 4 -20a/3 3 -4a/3 + 23a/3 2 -/? 4 + 7/3 3 - 10/3 2 )/((-l + 0) (-1 + a)) 

(B.33) 

T 4) i (a, /3) = - 2(10 a 4 /? 3 + 35 a 4 /3 + 10 a 3 /? 4 - 34 a 2 /? 4 + 35 a /3 4 - 34 a 4 /3 2 - 3 a /3 5 

- 3 a 5 /? + 2 a 5 /3 2 + 2 a 2 5 + 23 a 3 + 5 + a 5 - 24 a - 12 2 - 12 a 2 

+ 114 a 3 /3 2 - 86 a 3 /? + 114 a 2 3 - 152 a 2 /? 2 + 79 a 2 /? - 86 a 3 + 79 a /3 2 

+ 23 /3 3 - 58 a 3 3 - 12 a 4 - 12 /3 4 ) /((-l + a) 2 (-1 + /?) 2 ) 

(B.34) 

T 4fi (a,0) =220 2 + 22a 2 -5a 2 0-5a0 2 + 28a0-a 3 - 3 (B.35) 

55,2 (a, 0) = - {-a 2 + a 2 - a + a 2 - 2 ) (6 a + 6 - 6 a 2 - 6 2 - 30 a /3 - 36 a 2 2 
+ 32 a 2 + 32 a 2 + 3 + a 3 + 9 a 2 3 + 9 a 3 /? 2 - 9 a 3 - 9 a 3 /?)/((- 1 + 0) 
(-1 + a) (-/3 - a + a/3)) 

(B.36) 

5 5 ,i(a, /3) = - (-5 3 + 52 a 2 /3 2 - 35 a 2 3 + 26 a 3 + 26 a 3 /? - 19 a 2 + 6 a 2 /? 4 
- 35 a 3 /3 2 - 5 a 3 + 6 a 4 2 - 19 a 2 /? + 2 a 4 - 7 a 4 /? + 12 a 3 3 + 2 A 
-7a0 4 )/(-0-a + a0) 

(B.37) 

5 5) i(a, /?) =(-133 a 4 /? 3 - 89 a 4 /? - 133 a 3 /3 4 + 176 a 2 /3 4 - 89 a 4 + 176 a 4 2 + 18 a 3 /3 5 
+ 18 a 5 /3 3 + 17 a 5 + 17a 5 /?- 35 a 5 /? 2 - 35 a 2 5 -12 a 3 - 319 a 3 /? 2 
+ 136 a 3 - 319 a 2 /3 3 + 236 a 2 /? 2 - 60 a 2 /? + 136 a /3 3 - 60 a 2 - 12 /3 3 
+ 310a 3 /3 3 + 36a 4 /3 4 + 10a 4 + 10/3 4 )/(2 (-1 + 0) (-1 + a) (-/3-a + a/?)) 

(B.38) 

,§5,0 (a, 0) = - (5 a 4 3 + 6 a 4 /? + 5 a 3 /? 4 - 20 a 2 /? 4 + 6 a /3 4 - 40 a 3 3 + a 5 2 - 20 a 4 /3 2 
-2a/3 5 -2a 5 /3 + a 2 /3 5 + 28a 3 /? 2 + 28a 2 3 + 2 5 + 2 a 5 )/(2 a /3 
(-/?-a + a/3)) 

(B.39) 

D 5) i(a, /?) = - (-178 a 4 /? 3 + 2 a 4 - 178 a 3 /3 4 + 112 a 2 /3 4 + 2 a /3 4 + 112 a 4 2 

- 28 a 3 /? 5 - 28 a 5 /3 3 - 33 a /3 5 - 33 a 5 /? + 45 a 5 2 + 10 a 5 /3 4 + 10 a 4 /3 5 
+ 45a 2 5 + Ua 3 6 + Ua 6 3 + 13a 6 + 13a/3 6 - 27a 2 /3 6 -27a 6 /3 2 
+ 6/3 5 + 6 a 5 - 200a 3 2 + 24 a 3 - 200 a 2 /3 3 + 64a 2 /3 2 + 24a/3 3 

+ 370 a 3 /3 3 + 62 a 4 /3 4 - 8 a 4 - 8 /? 4 )/(2 (-1 + /?) (-1 + a)(-0-a + a 0)) 

(B.40) 



- 30 - 



T 5)2 (a, (3) =(18 a 5 /3 d - 24 a 5 /3 2 + 5 a 5 /? + a 5 + 36 a 4 /T - 84 a 4 /? 3 + 78 a 4 /3 2 - 25 a 4 (3 



a 4 + 18a 3 /3 5 + 98a 3 /3 3 



60 a 3 f3 2 + 24 a 3 (3 - 84 a 3 /3 4 - 24 a 2 (3 5 



+ 78 a 2 (3 A + 2a z /? z -60a z + 5 a /3° - 25 a /3 4 + 24 a /3 d 
+ /3 5 - /3 4 )/(2 a /3 (-! + «)(-! + /?)) 



.2 a 2 



.2 «3 



(B.41) 



T 5> i (a, /?) = - (12 a 5 /? 3 - 20 a 5 (3 2 + 9 a 5 /? - a 5 + 24 a 4 (3 4 - 88 a 4 /? 3 + 94 a 4 /? 2 
- 31 a 4 /? + a 4 + 12 a 3 (3 b + 178a 3 /? 3 - 124a 3 /3 2 + 24a 3 - 88a 3 /? 4 



- 20 a z (3° + 94 a 2 /3 4 + 46 a 2 /3 2 - 124 a 2 /? 3 + 9 a /? 5 - 31 a /3 4 



- 2 /3 5 

5 3 o5 



+ 24a (3 d - (3 b + (3 A )/(a /?(-! + a) (-1 + /?)) 



(B.42) 



T 5i0 (a, /?) = - (7 f3 A + 6 a 4 /3 2 + 32 a 3 /? + 6 a 2 /3 4 + 32 a /3 3 + 12 a 3 /3 3 - 14 a 4 /? + 7 a 4 
+ 66a 2 /3 2 - 70a 3 /? 2 - 14a/3 4 - 70a 2 /? 3 )/(2a/3) 

(B.43) 



C. Real-emission coefficients at NLO and NNLO 

Below we list the real-emission coefficients c[ n \a, (3 = ra) in Eq. (4.1) at NLO (n = 1) 
and NNLO (n = 2) of each of the five structure functions, i = 1 to 5. Note that these 
expressions include a singular (non-integrable) piece for r — > 0. According to the default 
subtraction prescription we use, Eq. (6.4), the plus prescription is defined with respect to 
r = (3 /a. The coefficients c\ n \a,(3 = ra) are therefore understood to be separated as 
in Eq. (6.6), where the the singular part that is put under the plus prescription as given 
in Eq. (6.7). The NLO coefficients are: 



c\ (a, ar) 



(2a 2 - 10a + 7) r 3 + (2a 2 - 4a - 7) r 2 - (10a + 7)r + 7 

4(r - l)r 

(r 4 + r 2 ) a 2 - 2r (r 3 + 2r 2 + 2r + l) a + 2 (r 2 - l) 2 



2(r - l) 2 r 



lnr 



(C.l) 



c 2 (a, ar) 



(-6a 3 + 20a 2 - 21a + 7) r 3 - (6a 3 - 56a 2 + 57a - 10) r 2 



a(r — l) 3 

(20a 2 - 57a + 38) r (10 - 21a)r + 7 
~b 7 7Tq ~b 



a(r — l) 3 



a(r — l) 3 r 



+ 2 



a 



- 3a 2 + 4a - 2) r 4 + (4a 3 - 21a 2 + 18a - 2) r 3 



a(r - l) 4 

(a 3 - 21a 2 + 34a - 14) r 3 + (-3a 2 + 18a - 14) r 2 + (4a - 2)r - 2 

a(r — l) 4 r 



lni 



(C.2) 



- 31 - 



Cg 1 ^ (a, ar) 



(4a 2 - 10a + 7) r 2 + (2 - 10a)r + 7 
2a(r — l)r 

(a 2 - 2a + 2) r 3 + (a 2 - 6a + 2) r 2 - 2(a - l)r + 2 
a(r — l) 2 r 



(C.3) 



lnr 



c\ (a, ar) 



+ 



22r 2 + 28r - a (r 3 + 5r 2 + 5r + l) + 22 
a(r-l) 3 

4 ((a - 2)r 3 + (a - 7)r 2 + (a - 7)r - 2) 
a(r-l) 4 



(C.4) 



lnr 



(-6a 2 + 14a - 7) r 3 - 2 (6a 2 - 35a + 16) r 2 - (6a 2 - 70a + 66) r 



c 5 (a, ar) 



2(7a - 16)r - 7 
2a(r — l) 3 r 



a 



2a(r - l) 3 
2 - 2a + 2) r 4 + (5a 2 - 20a + 6) r 3 



a(r - l) 4 



(C.5) 



(5a 2 - 40a + 28) r 2 + (a 2 - 20a + 28) r -2(a - 3)r + 2 

"I" 7 -i \ A l~ 



a(r-l) 4 

The NNLO coefficients in the large-/3o limit are: 



a(r — l) 4 r 



lnr 



c- (a, ar) 



-1 



s , / n ™ , x a(l — r)(l + r — ar) , , 
Ai(a, r) ln(r) + Bj(a, r) -i ^ ^ ln(l + r - ar) 



" [a(l - r)f m 
+ Tj(a, r) + 2a(l + r - ar) S^a, ar 
1 



Li 2 



r 

1 — ar 
1 + r — ar 



- Li 2 



r (1 — a) 
1 + r — ar 



+ — ln 2 (r) — ln(r) ln(l + r — ar 



+ ( ^ - 21n(a) -ln(r) ) c\ 1] (a,ar) 



(C.6) 



with 



Ti(a,r) =a(l-r) 



T if i(a, ar) + -T i)2 (o;,ar) 



where T i: j(a, (3) and Si^a, (3) in (C.6) are given in Appendix B and yi = [1, 2, 1, 2, 2]. The 
other functions Aj(a,r) and Bj(a,r) are: 

„2 



vi(a,r) = 



a 



(2a 4 - 4a 3 + 3a 2 ) r 6 + (2a 4 - 7a 3 + 6a 2 - 5a) r 5 



2r(ar — l) 2 

+ (a 4 - 11a 3 + 8a 2 + 5a + 2) r 4 + (l9a 2 - 8a - 6) r 3 
+ (-5a 2 - 12a + 4) r 2 + (6a + 2)r - 2 



(C.7) 



V2(a,r) = a 3 



-8a 3 + 20a 2 - 18a + 6) r 4 + (-20a 3 + 78a 2 - 84a + 26) r 3 



- (2a 3 - 36a 2 + 76a - 36) r 2 + (-2a 2 - 16a + 20) r + 8a - 



(C.8) 



- 32 - 



A 3 (a,r) 



a 



(3a 3 - 4a 2 + 3a) r 4 + (a 3 - 8a 2 + 6a - 2) r 3 



r(ar — 1) 
+ (a 2 + 6a - 4) r 2 - 4ar + 2 



(C.9) 



t(a,r) 



or 



— aV - (10a 4 - 21a 3 + 4a 2 ) r 6 - (l3a 4 - 64a 3 + 68a 2 - 6a) r 5 



(ar — l) 3 

+ (-6a 4 + 67a 3 - 138a 2 + 72a - 2) r 4 + (lOa 3 - 116a 2 + 136a - 26) r 3 
+ (8a 2 + 76a - 48) r 2 + (-20a - 16)r + 8 



(CIO) 



A 5 (a,r) 



a 



-Aa A + 6a 3 - 3a 2 ) r 7 + (-14a 4 + 47a 3 - 35a 2 + 8a) 



r (ar — l) 2 

+ (-11a 4 + 80a 3 - 132a 2 + 52a - 4) r 5 + (-a 4 + 29a 3 - 133a 2 + 134a - 26) r 4 
+ (-20a 2 + 84a - 42) r 3 + (5a 2 - 2a - 18) r 2 + (4 - 6a)r + 2 



(C.11) 



a 



Ma,r 



(-a 5 + 7a 4 - 11a 3 + 5a 2 ) 



2(a- l) 2 (ar-l) 2 

+ (-a 5 + 4a 4 - 10a 3 + 17a 2 - 9a) r 3 + (7a 4 - 10a 3 - 6a 2 + 3a + 4) r 2 (C.12) 
+ (-11a 3 + 17a 2 + 3a - 8) r + 5a 2 - 9a + 4 

h(a, r) = 2a 2 [ (3a 2 - 8a + 5) r 3 + (3a 2 - 26a + 13) r 2 + (13 - 8a)r + 5] (C.13) 



a,r) = 



(-2a 3 + 6a 2 - 5a) r 2 + (6a 2 - 8a + 4) r - 5a + 4 
(a — l)(ar — 1) 



(C.14) 



i(a,r) 



a 



(a-l) 3 (ar-l) 3 



a 



2a 5 + a 4 ) r 6 + (5a 6 - 38a 5 + 68a 4 - 37a 3 + 4a 2 ) r 5 



+ (5a 6 - 52a 5 + 171a 4 - 230a 3 + 104a 2 - 6a) r 4 

+ (a 6 - 38a 5 + 171a 4 - 306a 3 + 282a 2 - 100a + 2) r 3 

- (2a 5 - 68a 4 + 230a 3 - 282a 2 + 160a - 34) r 2 

+ (a 4 - 37a 3 + 104a 2 - 100a + 34) r + 4a 2 - 6a + 2 



(C.15) 



- 33 - 



a 



(a - l) 2 (ar - 1) 



a - l)V(3a - 5)r J 



+ a (6a 4 - 49a 3 + 87a 2 - 55a + 12) r 4 

+ (3a 5 - 49a 4 + 160a 3 - 180a 2 + 68a - 6) r 3 

+ (-11a 4 + 87a 3 - 180a 2 + 140a - 30) r 2 

+ (13a 3 - 55a 2 + 68a - 30) r - 5a 2 + 12a - 6 



(C.16) 



D. Virtual coefficients at NLO and NNLO 

The coefficients v^' 2 \a) entering (7.21) for the five structure functions, i = 1 to 5 are 
listed below. These coefficients are computed by expanding Eq. (7.20). Let us recall that 
these expressions correspond to the infrared-subtraction procedure detailed in Sec. 7, where 
£[Vfo](a:,ii)| si is defined according to (6.12), where the plus prescription is defined with 
respect to r = (3 /a, as in (6.4). 

Let us also recall that according to Eq. (7.20) there is a simple all-order relation 
between the virtual corrections for the structure functions i = 3 and i = 1, namely 



in) 



(a) 



a 



a) . 



The NLO and NNLO coefficients for the other structure functions are: 



(D.l) 



«J 1} (a) 



2Li 2 (l-a) 



(-3 + 2a)ln(a) 2vr 2 



-1 + a 



(D.2) 



(2), s a 
v \ (") = g 



~ 21n(l - a)) ln 2 (a) + ( §^ + 2^ 2 ) ln(a) 



+ 



a — 1 
16 - 19a) Li 2 (l - a) 
3(a - 1) 



6(a - 1) 



2Li 3 (l-a)-4Li 3 (a) + 2C 3 



79 vr 2 71 



36 



12 



(D.3) 



v£\a) = -41n(a) - 4Li 2 (l - a) - - 5 



47T 2 

3 



(D.4) 



4 2) (a) = (4-41n(l-a))ln 2 (a) + I 4tt 



2 (7a - 4) 
3(a - 1) 



ln(a) 



38 797r 2 71 

- yLi 2 (l - a) - 4Li 3 (l - a) - 8Li 3 (a) + 4( 3 - — - — 



(D.5) 



^ 1} (a) 



2(-l + 2a)ln(a) 
(-1 + a) 2 



+ 



-1 + a 



(D.6) 



- 34 - 



(2) . 2(-l + 2a)Li 2 (l-a) 2 (-1 + 2 a) ln(a) 5 

U 4 («) = r-T, X2 1 ( 1 i \2 

(— 1 + a) z (— 1 + a) z 

1 (-25 + 38a)ln(a) ^ 19 



(D.7) 



3 (-1 + a) 2 3(-l + a) 



(i). . nT . „ , (2a 2 -2a + l)ln(a) 2tt 2 5q-7 q , 



)/ x / , / 2(a-l)a + l\, , /14a 2 + 16a -17 9 \ , . . 

22 + a(19a-44) T . ^ , oT .. i . , T . . , 109 - 71a 79^ 2 
+ 3(a l _ 1)2 L l2 (l - a) + 2L 13 (1 - a) + 4L l3 (a) + + — 

(D.9) 



E. Polynomials entering the triple differential width 

Below we list the polynomials in r = (3 /a, a and x entering the triple differential width 
summarized in Sec. 8. The polynomials entering the NLO coefficient through Eq. (8.10) 
are: 

Qi(a, r, x) = a 2 (1 - a) (a 2 - 2 a + 2) r 6 + a(l - a) x 

(-7 a 2 + 2xa 2 + 5a-4xa-6 + 4x)r 5 - (-10 x a + 42 a 3 - 2x 2 + x 2 a 3 
+ 29 x a 2 - 4 + 4 ax 2 + 4xa 4 - 3x 2 a 2 + 4a 5 - 20a 4 + 6x + 3a - 28 x a 3 

- 25 a 2 )r 4 - (-64 x a + 4 - 98 a 2 + 4x a 4 - 21x 2 a 2 + 46 a - 40 x a 3 

- 4x- 20a 4 + 95xa 2 + 18ax 2 + 68a 3 + 4x 2 a 3 )r 3 - (75xa 2 - 80xa (E.l) 
+ 34 ax 2 + 32 a 3 + 52 x + a 5 - 21 x 2 a 2 + x 2 a 3 - 22 x 2 + 52a + 2xa 4 - 7a 4 

- 48 a 2 - 30 - 28 x a 3 )r 2 + (- 16 a - 18 a x 2 + 2 x 2 + 3 x 2 a 2 - 10 - 39 x a 2 

- 20 a 3 + 43 a 2 + 8 x + 6 x a 3 + 44 x a + 3 a 4 )r - 2 a 3 - (6 x - 13) a 2 

- (27 + 4x 2 - 26x)a + 16+ 10x 2 - 26x 



a, r, x) = 2a 2 r 2 - 10 a r 2 + 7r 2 + 2a 2 r - 4a r - 7r - 10a (E.2) 



h(a, r, x) = (6a 3 - 20a 2 + 21a -7) r 3 + (-17 - 56a 2 + 57a + 6a 3 ) r 2 
+ (-20 a 2 - 17 + 57 a) r + 21 a - 31 



(E.3) 
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The polynomials entering the NNLO coefficient through Eq. (8.11) are: 

Pi (a, r, x) = a 3 (1 - a) (3 a 2 - 12 a + 10) r 6 

- a 2 (38 + 3 a 4 - 20 x - 38 x a 2 + 48 x a - 22 a 3 - 83 a + 63 a 2 + 10 x a 3 ) r 5 

- a(26 a x 2 - 46 + 3 a 5 + 6 x 2 a 3 + 4 x a 4 - 10 x 2 - 58 a 4 - 140 x a - 82 x a 3 

- 160 a 2 - 22 x 2 a 2 + 170 x a 2 + 155 a 3 + 110 a + 48 x)r 4 - (-334 a 3 - 28 x 
+ 155a 4 + 10x 2 - 22a 5 + 136xa - 68a - 64x 2 a 3 + 242a 2 - 52ax 2 

+ 284 x a 3 - 82 x a 4 - 320 x a 2 + 6 x 2 a 4 + 100 x 2 a 2 + 3 a 6 + 10 x a 5 + 18)r 3 
+ (15 a 5 + 22 x 2 a 3 - 63 a 4 + 320 x a 2 - 26 x 2 + 44 x + 104 a x 2 + 38 x a 4 

- 100 x 2 a 2 + 160 a 3 - 232 x a + 144 a - 242 a 2 - 18 - 170 x a 3 )r 2 - (136 x a 

- 140 x a 2 + 26 x 2 + 26 x 2 a 2 - 44 x + 18 + 22 a 4 + 110 a 2 - 52 a x 2 - 68 a 

- 83 a 3 + 48 x a 3 )r + 10 a 3 + 2 (- 19 + 10 x) a 2 + 2 (23 - 24 x + 5 x 2 ) a - 18 
+ 28x - 10x 2 

(E.4) 

P 2 (a, r, x) = -2 a 3 (1 - a) (34 a 2 - 74 a + 59) r 7 

- a 2 (1 - a) (24 a 3 - 424 a 2 + 184 x a 2 + 557 a - 344 x a - 460 + 236 x) r 6 

+ 2a(289 x + 650 a - 283 - 1204 a 2 + 1229 x a 2 + 1291 a 3 - 822 x a + 52 x 2 a 3 

- 804 x a 3 - 539 a 4 + 157 a x 2 + 76 x a 4 - 59 x 2 - 150 x 2 a 2 + 85 a 5 )r 5 

+ (3320 a 4 - 1208 x 2 a 3 + 224 - 667 a + 200 x 2 a 4 - 5680 a 3 - 342 x - 4674 x a 2 

+ 5352 x a 3 + 118 x 2 - 554 a x 2 + 3451 a 2 + 1416 x a - 1732 x a 4 + 24 a 6 

+ 260 x a 5 + 1380 x 2 a 2 - 672 a 5 )r 4 + 2(2592 a 2 + 56 + 2 x a 5 - 127 a 5 

+ 1936 x a 3 + 727 a 4 - 470 x a 4 + 1120 x 2 a 2 - 3424 x a 2 - 2 x 2 a 4 - 2105 a 3 

- 176 x + 7 a 6 + 120 x 2 - 1150 a + 1910 x a - 820 a x 2 - 310 x 2 a 3 )r 3 + (16 a 5 

- 1076 x + 76 x 2 a 3 - 1524 a x 2 + 560 x 2 + 660 x 2 a 2 - 1491 a 3 + 104 x a 4 

- 2178 a + 944 x a 3 + 2834 a 2 - 3424 x a 2 + 303 a 4 + 516 + 3744 x a)r 2 
-2(412x + 37a 4 + 79ax 2 + 83x 2 a 2 - 218- 194x 2 - 16xa 2 + 123xa 3 -493xa 
+ 401 a - 184 a 2 - 36 a 3 )r + 44 a 3 + (-145 + 138 x) a 2 

+ (-296 x + 94 x 2 + 165) a - 64 - 82 x 2 + 146 x 

(E.5) 

P 3 (a, r, x) = a (1 - a) (50 a 2 - 142 a + 85) r 3 - (50 a 4 + 219 a 2 - 142 a + 85 - 200 a 3 ) r 2 
+ (-304 a 2 + 192 a 3 + 39 a + 85) r - 142 (- 1 + a) a 

(E.6) 

P 4 (a, r, x) = (1 - a) (150 a 2 - 218 a + 85) r 3 - (-1016 a 2 + 1035 a - 287 + 150 a 3 ) r 2 
+ (368 a 2 - 1035 a + 395) r + 457 - 303 a 

(E.7) 



- 36 - 



P 5 (a, r, x) = -(85 + 88a 2 - 142 a) r + 142 a - 123 



(E.8) 
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